Topologically protected Landau levels in bilayer graphene in finite electric fields 
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The zero-energy Landau level of bilayer graphene is shown to be anomalously sharp (delta-function 
like) against bond disorder as long as the disorder is correlated over a few lattice constants. The 
robustness of the zero-mode anomaly can be attributed to the preserved chiral symmetry. Unex- 
pectedly, even when we apply a finite potential difference (i.e., an electric field) between the top and 
the bottom layers, the valley-split n = Landau levels remain anomalously sharp although they are 
now shifted away from the zero energy, while the n = 1 Landau levels exhibit the usual behavior. 



Introduction — The existence of the zero-energy Lan- 
dau level is a most fundamental property of the electronic 
states in graphene in magnetic fields, which hallmarks the 
unconventional quantum Hall effect observed in mono- 
layer grapheneJ^ Specifically, the zero-energy Landau 
level of mono-layer graphene shows an anomalous robust- 
ness against the disorder induced by ripples, an intrinsic 
disorder in graphene, which leads to an unconventional 
criticality of Hall transition at zero energy^ii For the ro- 
bustness of zero modes, the chiral symmetry^ defined in 
terms of the chiral operator T that anti-commutes with 
the Hamiltonian H, {T, H} = with T 2 = 1, is an essen- 
tial ingredient. 6 For monolayer graphene, we have ver- 
tical Dirac cones at K and K' points in the Brillouin 
zone, and the effective Hamiltonian has the chiral sym- 
metry. In such a system, it has been demonstrated by 
the present authors that the zero-energy (n = 0) Landau 
level is robust against the disorder that respects the chi- 
ral symmetry as long as the disorder is correlated over a 
few lattice constants^ Experimentally, the n = Lan- 
dau level narrower than the other n ^= Landau levels 
is reported for monolayer graphene,— which is consistent 
with the present robustness specific to the zero-energy 
(n = 0) Landau level. 

The notion of the chiral symmetry is so universal 
that it has further been showniS that the chiral sym- 
metry, usually considered for the vertical Dirac cones, 
can be generalized to accommodate tilted Dirac cones, 
such as those encountered in certain organic metalsJ^— 
The generalized chiral symmetry protects the zero-energy 
Landau level as far as the Hamiltonian as a differential 
operator is elliptic, where we can even extend the argu- 
ment of Aharonov and Casher for counting the number of 
zero modes in the presence of disorder^ The existence 
of the generalized chiral symmetry can indeed be trans- 
lated to a condition that the index theorem-^ holds for 
generic tilted Dirac cones. The chiral symmetry is there- 
fore directly related to the robustness of zero modes for 
the massless Dirac fermions. 

Now, in the physics of graphene, the case of bilayer 
graphene is an interesting test bench for examining vari- 
ous graphene properties. Specifically, McCann and Falko 



have shown that there exist four-fold degenerated (per 
spin) zero-energy Landau levels, which lead to a quan- 
tum Hall effect characteristic to bilayer graphenej** The 
degeneracy comes from the valley (K and K') degrees of 
freedom and two (n = and n = 1) Landau indices. Al- 
though the robustness of these zero mode Landau levels is 
also predicted as a consequence of the index theorem ) 17 i 18 
it is not clear whether there is also a direct relationship 
between the chiral symmetry and the anomalous robust- 
ness of zero modes, since the parabolic band dispersion 
in the bilayer graphene, as opposed to the linear one in 
the monolayer graphene, might well invalidate the argu- 
ments. For instance, the robustness of the n = level 
and that of n = 1 level can naively be different, since 
they have different structure for the wave functions. 

The difference becomes even greater when we apply an 
electric field perpendicular to the graphene sheet, which 
introduces an energy gap. An opening of the energy-gap 
in bilayer systems is important in an applicational con- 
text as welli2i~— Thus a further interest is to see what 
happens to the zero-energy Landau level when the en- 
ergy gap is introduced. Experimentally, it is desirable to 
clarify quantitatively the robustness of the zero modes 
in bilayer graphene, since experimental results in high 
mobility samples are now availablej 19 ' 20 

The purpose of the present paper is to explore these 
very questions, for which we have performed numerical 
studies based on the lattice model. We shall show that 
both the n = and the n = 1 Landau levels in bilayer 
graphene are robust against bond disorders as long as 
they are correlated over a few lattice constants. We anal- 
yse the result in terms of the chiral symmetry for the bi- 
layer system. Unexpectedly, it is further found that, even 
in the presence of a potential difference (an electric field) 
between the top and the bottom layers, the n = Lan- 
dau levels remain robust although they are shifted away 
from zero energy, while the n — 1 Landau levels exhibit 
the usual behavior. This phenomenon is also discussed 
in terms of the effective theory at K and K' points. 

Lattice model — In order to investigate the robustness 
of the zero modes against disorder in bilayer graphene, 
we adopt the following tight-binding lattice model with 
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FIG. 1: (Color Online) For a bilayer graphene we show the 
A-B stacking with the interlayer coupling 71 (left panel), and 
the band dispersions (right). K — (2tv/S\/3, 2ix/i)a~ 1 and 
K' = (-27r/3V3,27r/3)a _1 are the corners of the Brillouin 
zone. 



the Bernal (A-B) stacking ji 6 . We assume that each layer 
can be described by the simple honeycomb lattice, while 
the interlayer coupling, 71, connects a site (Bi) on the 
B sublattice of the bottom layer and a site (A2) on 
the A sublattice of the top layer just above Bi(FigJTJ 
left). This simplest model accounts for the parabolic 
dispersion with zero gap at K and K' points of bilayer 
graphene (Fig[TJ right). For the randomness, we con- 
sider a bond disorder that is spatially correlated. This 
is described by a random component, 5t(r), for the hop- 
ping amplitude in each layers as t(r) = t + St(r), that is 
gaussian-distributed with a variance a and is correlated 
in space with a correlation length rj as (St(r)St(r')) = 
(St 2 ) exp(— \r — r'| 2 /4?7 2 ). It is to be noted that in this 
tight-binding lattice model the chiral symmetry is exactly 
preserved even in the presence of disordered components 
in hopping amplitudes t and Jir^ which is due to the 
bipartite structure of the lattice. Spin degrees of freedom 
are suppressed for simplicity. 

A randomness is expected to be induced in the hopping 
amplitudes as a consequence of ripples 2 ^ in a monolayer 
graphene. In the case of bilayer, the disorder should be 
correlated between the two layers if the two layers have 
a common ripple. In present paper, however, we also 
examine the case where the disorder in two layers are 
uncorrelated to clarify the generality of the topological 
protection of the Landau levels in bilayer graphene. 

The effect of the magnetic field is taken into account 
by the Peierls substitution t —> te~ 27Tl8 ( r \ such that the 
summation of the phases along a loop is equal to the 
magnetic flux enclosed by the loop in units of the flux 
quantum </>o = (h/e). The nearest- neighbor distance of 
the honeycomb lattice is denoted by a, while the external 
uniform magnetic flux enclosed by the hexagon of the 
honeycomb lattice by 4>. 

Zero-mode Landau level — Let us first discuss the ro- 
bustness of the zero-energy Landau levels. The density 
of states p(E) = —(lmGu(E + i£)/ir)i with Gu(E + ie) — 
(i\(E — H + ze) _1 |i) is evaluated by the Green function 
method^ Figure [5] displays the density of states for the 
case where the disorders in the two layers are perfectly 




FIG. 2: (Color Online) Density of states of bilayer graphene 
in a magnetic field plotted for varied correlation length, 77, 
of the bond disorder. The result is for a system with 10 6 
sites, a magnetic field 4>/4>o = 1/50, the degree of disorder 
a/t — 0.115, inter-layer transfer 71/t = 0.2, and e/t = 0.0006. 
Energy dispersion in zero magnetic field around the K point 
is also shown with the same energy scale. Inset: Schematic 
figure of correlated ripple. 

correlated. We find that the zero energy Landau level 
becomes anomalously sharp as soon as the spatial corre- 
lation length 77 exceeds a few nearest-neighbor distances 
a, which is the same behavior as in the case of the mono- 
layer graphene. This means that both of the n = and 
the n = 1 Landau levels at zero energy remain delta- 
function like in the presence of finite-range bond disor- 
der, which confirms the prediction based on the effective 
Hamiltonian at low energies^ We also examine the case 
where the disorder in two layers is uncorrelated to find 
that the density of states coincides with those in the cor- 
related case, which reveals an insensitivity to the disorder 
correlation between two layers (Figj3]). 

Electric field effect — In a bilayer graphene we can in- 
troduce an energy gap by applying an electric field per- 
pendicular to the graphene sheet (FigHJ inset). 23 While 
the electric field obviously breaks the inversion symme- 
try of the system, we should note that it also breaks the 
chiral symmetry for the effective Hamiltonian around K 
and K'. To examine what happens to the robustness of 
the Landau levels around E = in such a case, we evalu- 
ate the density of states when the potential difference 2A 
is introduced between the two layers, where the four-fold 
degenerated zero-energy Landau levels split into fourJ£ 
The present numerical result, displayed in Fig 01 shows 
that, among the split four Landau levels, the n — Lan- 
dau levels that are located at energies E = ±A again 
become anomalously sharp as soon as the disorder is cor- 
related over few lattice constants, while the n = 1 Lan- 
dau levels are broadened by the disorder despite the fact 
that they are located closer to E — 0. An interesting ob- 
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FIG. 3: (Color Online) Density of states of bilayer graphene 
for the case of correlated disorder between two layers (dotted 
curves) and that for uncorrelated disorder (solid curve) are 
shown with an offset 1.0 for uncorrelated disorder. The other 
parameters are the same as in Fig[5] 




FIG. 4: (Color Online) Density of states of bilayer graphene 
in the presence of a potential difference A/t = 0.1. The other 
parameters are the same as in Fig. Q] Energy dispersion in 
zero magnetic field around the K point is also shown in the 
same energy scale. 



servation is that the energies of these anomalously sharp 
Landau levels can be tuned by the electric field. We can 
also note that the n = 1 Landau levels, while not as sharp 
as the n — Landau levels, are significantly sharper than 
higher Landau levels. We have also confirmed that the 
results are insensitive to whether the disorder is corre- 
lated between the two layers or not (FigJS]). 

Effective theory — The effective Hamiltonian, act- 
ing on the envelop functions , ^ , W% i *b 2 ) at K " 
point, for bilayer graphene in a perpendicular electric 
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FIG. 5: (Color Online) Four Landau levels split by the elec- 
tric field around zero energy of bilayer graphene when the 
disorder correlation length in each layer is r//a — 2. Here 
A/t — 0.1 and the other parameters are the same as in Fig. 
[2] The results (solid curve) for the case where the disorder 
is uncorrelated between two layers is plotted on those(dotted 
curve) for the case of the same disorder in two layers. Two 
curves almost coincide, demonstrating clearly the insensitiv- 
ity to whether the disorder is correlated between the layers or 
not. A small width of the n — levels is an artifact of a finite 
imaginary energy e in the Green function. Inset: The energies 
of the n = 1 Landau levels in the presence of disorder are plot- 
ted by solid circles as a function of the interlayer coupling 71 . 
The solid(dotted) curve represents the perturbational results 
(Eq. CO)) for K(K') . 
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where ir? = p x £ — ip y i with pi = — ihd + eAi{e > 0), and 
vf the Fermi velocity of a monolayer graphene. Here 
£(= 1,2) labels the two layers, and the gauge field An 
represents the effective gauge field in each layer, which 
includes the random gauge field induced by the random 
hopping as well as the contribution by the external (uni- 
form) magnetic field. The operator ixg satisfies a commu- 
tation relation, [7T^,7r!] = 2heBi, with Bi — (V x Af) z . 
Note that the gauge fields A\ and A2, and consequently 
the effective magnetic fields B\ and B2 at K point, can be 
different when the random bonds are different between 
the two layers. 

An important point is that the random component in 
the gauge field induced by the bond disorder gives rise 
to effective magnetic fields that have opposite signs be- 
tween K and K' points^ 3 Still, the effective Hamiltonian 
is chiral-symmetric, since it satisfies THT = — H as long 
as A = 0, where the chiral operator L is given, for a 
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bilayer, in terms of a Pauli matrix a z as 
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The zero modes for the valley K with no electric field 
(A = 0) are, in analogy with the case without disorder, 25 
given by 
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with ttiV'o = 0' ""l^i — V'o 2 '' and ^ipo = 0. Note 
that these zero modes are also eigenstates of the chiral 
operator T. 

Following the argument by Aharonov and Casher— , 
we adopt the Coulomb gauge d%A x + dyA y = 0, and ex- 
press the gauge field as A = (—d y (p, d x <f>). The operator 
7r is then expressed as m = —2ih[d z * + 2n(d z *<pe)/(j>o\ 

with z = (x + iy)/2. The solution to Treipjp — with 

I = 1,2 is given by % 1 ' = fe(z)exp(—2ir<pi/<f>o), where 
fl(z) a polynomial in zJ^ In general, the solution ipi 
takes the form ip\ — f2(z)exp(—2iripi/(j)o)F(z,z*), with 
d z *F(z,z*) — exp(2Tr(ipi — ip 2 )/<l>o). If the bond disor- 
ders in two layers are the same (ifi = if 2)1 the func- 
tion F(z,z*) is reduced to z*J£ It is straightforward to 
apply these arguments to the effective Hamiltonian for 
the valley K' ) 16 ' 25 i 26 which implies that the four- fold de- 
generated zero-energy Landau levels exist irrespective of 
presence or absence of the disorder in gauge fields. 

When the electric field is switched on (A =^ 0), the 
above chiral symmetry is broken. We can still show, how- 
ever, that the state ^n=o rem ains to be an exact eigen- 
state of the Hamiltonian with the eigenvalue e^" = —A, 
so that the broadening due to disorder is absent as in the 
case of A = is demonstrated. The state *„ = i, on the 
other hand, is not an exact eigenstate for A / 0. It is 
therefore natural to expect that the broadening occurs 
for the Landau levels corresponding to *^ = i, as is actu- 
ally seen in our numerical results (Figs. 0]and[S]). We 
can also note that even for such states, the broadening 
itself is likely to be significantly smaller than those for 
higher Landau levels, which comes from the anomalous 
character of the unperturbed Landau level ^^ = i (Fig. 

n. 

The eigenvalue Ef- for the state *&n=i can be esti- 
mated, in the absence of disorder, from the perturbation 
with respect to A, which gives 
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where B = B\ = B2 denotes the uniform external mag- 
netic field. Relations to the tight-binding parameters 



are given by v F = (3/2)at/h and cj) = (3\/3/2)Ba 2 . 
The energy Ef- is then estimated as Ef^ ~ 0.7 A for 
4> = (l/5O)0o an d Ji/t = 0.2, which accurately agrees 
with the present numerical result (FigjS]). The same ar- 
gument for the effective Hamiltonian at K' point leads 



to the Landau levels at E^ = —E^ = A and Ef 



—E^ ~ — 0.7A. We also examine the positions of the 
n = 1 Landau levels for various values of the interlayer 
coupling 71 to confirm that their positions are in good 
agreement with the perturbational result (Eq.(TT])) for the 
range < 71 < t (FigO inset). 

The peak heights of these valley-split Landau levels de- 
pend on the effective magnetic field for each valley. Since 
the effective fields induced by ripples in K and K' points 
have the opposite sign with the same magnitude^ the 
degeneracies of the Landau levels can be different for K 
and K', although their sum should be a constant^ In 
our numerical results, however, no significant difference 
in the peak heights for the valley-split Landau levels is 
seen (FigsH] and [5]). This can be attributed to the fact 
that the present density of states is an average over the 
sample. Our sample-size is much larger than the corre- 
lation length of bond disorder and the periodic bound- 
ary condition is adopted along the strip geometry^ It is 
therefore expected that the local fluctuation of the mag- 
netic field due to the bond disorder (ripples) is canceled. 

Conclusions — We have demonstrated, both numeri- 
cally and analytically, that the zero-energy Landau lev- 
els of bilayer graphene become anomalously sharp when 
the bond disorder is correlated over a few lattice con- 
stants. The anomaly is shown to be insensitive to the 
disorder correlation between the top and the bottom lay- 
ers, which suggests a relevance of the chiral symmetry 
to the present anomaly as in the case of the monolayer 
graphene. Another new finding is that the anomaly at 
the n = Landau level persists even in the case where 
the chiral symmetry for each valley is broken by the po- 
tential difference between two layers. The splitting of the 
pair of sharp Landau levels is controlled by the electric 
field, and their peak heights reflect the effective magnetic 
field strength at each valley. The anomalous sharpness 
of these levels found here may help to detect experimen- 
tally the local fluctuation of the effective magnetic field 
arising from ripples as unbalanced peak heights of these 
Landau levels in the local density of states. 
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